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Abstract 

Many zero temperature transitions, involving the deviation in the value of a 
U(l) conserved charge from a quantized value, are described by the dilute Bose 
gas quantum critical point. On such transitions, we study the consequences of 
perturbations which break the symmetry down to Zpj in d spatial dimensions. 
For the case d = 1, N = 2, we obtain exact, finite temperature, multicritical 
crossover functions by a mapping to an integrable lattice model. 
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The zero temperature (T), quantum phase transition in a dilute Bose gas has recently 
attracted some interest [|l]-0] because of its importance in a variety of different physical 
situations. For bosons with repulsive interactions in a chemical potential, /i, the quantum 
critical point is at \x = 0, where the T = density of bosons has a non-analytic dependence on 
/i. This quantum critical point controls the finite T quantum-classical crossovers in the dilute 
Bose gas in addition it is the critical theory for (?) the Mott-insulator to superffuid 

transition in a lattice boson system at a generic /i ||, (ii) the onset of uniform magnetization 
in a gapped quantum antiferromagnetic in a magnetic field ||, (m) the deviation from a 
saturated polarization in a quantum ferromagnet ||, and possibly other physical systems ||. 
A common feature of all these systems is that the Hamiltonian always has at least a global 
U(l) symmetry, and the transition involves a deviation in the expectation value of the U(l) 
conserved charge from a quantized (possibly zero) value. 

In this paper we will study the consequences of perturbations which break the U(l) sym- 
metry down to Zn (the cyclic group of N elements), and compute associated exponents and 
crossover functions. Such an analysis is useful in the experimental spin systems noted above 
(especially (ii) where crystal-induced spin anisotropy destroys the U(l) symmetry. Our 
most detailed results will be for N = 2 and spatial dimension d — 1, which is also the most 
interesting and non-trivial case: we shall obtain explicit results for multicritical crossovers 
as a function of /i, T, and the strength of the [/(l)-breaking perturbation. These are the 
first nontrivial, exact results for universal finite T crossovers near a quantum multicritical 
point in any system. 

We begin by reviewing the U(l) symmetric Bose gas theory. The theory is described by 
the action 

s = / A jT + ^IVAP - + , (i) 

where f b is a complex scalar field, m is the boson mass and we use h — ks — 1- A 
T = renormalization group analysis shows that the interaction u is irrelevant for d > 2 . 
Near d = 2, we define the dimensionless bare coupling go by w = (2m)K £ go/ Sd where k is 

2 



a renormalization momentum scale, Sd = 2n d / 2 / \(27i) d T(d/2)) is a phase space factor, and 
e — 2 — d. The renormalized theory is defined by a single coupling constant renormalization 
of go to g, and no other renormalizations are necessary. This renormalization leads to the 
/3-function 



where Ad = 1/2 + 0(e) is a known constant independent of g, with the value the 0(e) terms 
depending upon the precise renormalization condition. The result (0) is valid to all orders in 
g. The quantum critical point of So is at \i = 0, and has dynamic exponent z — 2, correlation 
length exponent z/ = 1/2 for all d ||. At T = and /i small, the boson density n = (|\1/b| 2 ) 
obeys [|| n = A6(fi)2mfi for d > 2 (A is a non- universal, cut-off dependent, constant), while 
for d < 2 we have n = C ( ^6 l (/i)(2m/i) d//2 , with Q a universal number which has non-trivial 
contributions at each order in e || due to the fixed-point interaction g = g* = ej Ad. 

We now break the symmetry down to Z^. The most relevant perturbation which ac- 
complishes this is 



where the coupling can be taken to be real, without loss of generality. For d > 2, the 
scaling dimension of Xn is simply its canonical dimension, d + 2 — Nd/2; so A2 is relevant for 
all d > 2 (dim(A2) = 2), A3 is relevant for 2 < d < 4, while Ajv >3 are irrelevant for all d > 2. 
For d < 2, we need to consider the renormalization of insertions; these are non-trivial 
because Sn creates or annihilates iV bosons, and these pre-existing bosons can then interact. 
The two-loop computation of this renormalization yields the scaling dimensions 



dim(Ajv) =i-N- (N 2 /2 - N + l)e + (N/2)(N - 1){N - 2) ln(4/3) £ 2 + 0(e 3 ) (4) 



Notice that the 0(e 2 ) term vanishes for N = 2. This is also true for all subsequent terms, as 
the only renormalization of ^ 2 B comes from a single series of ladder diagrams, and we have 
dim(A2) = 2 — e to all orders in e; there is no similar simplification for N > 3. Later, we 



dg/dhin = -eg + A d g 2 , 



(2) 




(3) 
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will verify the value of dim(A2) in d = 1 by an entirely different method. So A2 is relevant 
in d = 1. Evaluation of the series (^) at e — 1 predicts that A3<n<5 are irrelevant in d = 1, 
while the Xn>g are again relevant. This result for N > 6 is surely an artifact of the poorer 
accuracy of the series at larger N, and it is likely that all Xn>3 are irrelevant in d = 1. 

The A 2 perturbation is expected to drive the quantum phase transition in iS +iS 2 (accessed 
by the control parameter //) into the universality class of the transverse-field Ising model. 
This latter class has z = 1 and upper critical dimension d = 3. We will not discuss the 
details of the crossover between the U(l) symmetric, dilute Bose gas fixed point and the 
transverse-field Ising fixed point for general d here: we shall confine our attention in the 
remainder of the paper to the d = 1 case where both fixed points are below their respective 
upper critical dimensions. 

We now present a detailed analysis for d = 1, N = 2. It has been argued that the d = 1 
critical properties of the action Sq are identical to those of the \i = critical point in a 
dilute, spinless, non-interacting, Fermi gas [|J. The Bose field, ^b, is related to the Fermi 
field, typ, by a continuum Jordan- Wigner transformation: 

= exp (-in J dx'^ F {x')^ F {x') s j^ F {x) (5) 

(we are momentarily interpreting ^b, as operators). We can use this mapping to deduce 
the mapping of <S 2 into the fermionic theory. We point split ^f%(x) into ^ b{%)^ b{% + a), 
rewrite in terms of ^ p using @, and expand in powers of a. Retaining only low order terms, 
we obtain the following fermionic form for Sq + S2 in d — 1: 

where A 2 oc A 2 , and fpa complex Grassmanian field. The action Sp is quadratic in fermionic 
fields, and correlations of the \&p can be computed. In particular, it can be shown that all 
other terms derivable from Sq + 5 2 are irrelevant at the fi = 0, A 2 = multicritical point. 
So we expect the universal, multicritical crossover functions to emerge from an analysis of 
Sp. Note also, from simple power-counting, that dim(A 2 ) = 1, which agrees with our earlier 
result for dim(A 2 ), obtained by expanding to all orders in e. 
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The T = phase diagram of Sp is shown in Fig. [|. The z = 2, dilute Bose gas critical 
point is M, and it is the point of intersection of three second order phase transition lines; 
there is a gap to all excitations everywhere, except at M and along these three lines, (i) The 
line along A2 = has U(l) symmetry, and describes the d — 1 Bose gas with quasi-long 
range order. At sufficiently long scales and for /x > 0, this line is described by its own 
critical theory which has z = 1, is conformally invariant and has central charge c = 1. 
There is an operator (l^p^^/^l 2 ) which is marginal along this line, and is responsible for 
the continuously varying exponents at c = 1; however this operator is irrelevant at the 
critical end-point M [|J and can be ignored while computing the multicritical crossovers of 
M. Under this condition, the c = 1 theory can be written as two copies of the c = 1/2 Ising 
field theory, (ii) The lines \i = 0, A2 > and \i = 0, A2 < are also conformally invariant 
at long scales, and are then described by a single z = 1, c = 1/2 Ising field theory. The 
non-zero expectation values for (^b) (Fig [1]) appear only at T = 0, and we always have 
(*&b) = f° r T > 0; this will become clear from our computations below. 

We now wish to describe the finite T crossovers in the vicinity of M. We will study the 
two-point correlators of ^b- It is useful to define ^ x = ^ b + ^b an d ^>y = —i{^B — *b)- 
Then, elementary considerations show that, for A2 real, (^>x^y) — and {^x^x)\ — 
(^y^y) t . So it is sufficient to compute i^x^x) for both signs of A2. We will describe 
the long distance behavior of its equal time correlation, where we expect 

lim (#x(x, r)* x (0, r)) ~ Ae~ lx ^ for T > 0. (7) 

\x\— *oo 

The correlation length, £, and the amplitude, A, obey the multicritical scaling forms 

r 1 = (2mT) 1/2 F(x,y) , A = 2 (2mT) 1/2 G(x,y) (8) 
where F and G are fully universal scaling functions of the dimensionless variables 

x = fi/T ; y = (2m) 1 / 2 A 2 /^ 1/2 (9) 

The powers of T in (^|,^|) follow from the exponents and scaling dimensions at M: z = 2, 
v = 1/2, dim(/i) = 2, dim(A2) = 1, and dim(\E , s) = 1 (the mass, m, is not to be interpreted 
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as a scaling variable; it converts between the engineering dimensions of space and time, and 
is analogous to the velocity of light in a Lorentz invariant theory). 

We will now provide an exact, closed-form, computation of the functions F and G. Our 
strategy is to perform the computation in an integrable lattice model with a multicritical 
point in the universality class of M. It turns out that the well-known Lieb-Schultz-Mattis ]7J 
spin chain has two such critical points. This spin chain is described by the Hamiltonian 

H = -J2{J((1 + 7KX + i + (1 - 7KV + i) + A*?} (10) 

i 

where a x,y ' z are Pauli matrices on the sites i of an infinite chain, J > 0. The Jordan- 
Wigner transformation maps H into a model of free, spinless, lattice fermions, and its phase 
diagram can then be computed exactly 0H1]; the result is shown in Fig. ||. The points Mi, 
M.2 are both in the universality class of M, and the continuum limit of the Jordan- Wigner 
transform of H yields precisely Sf] near Mi we find m = Ti 2 /{AJa 2 ), A 2 = 4:Jja, li = 2A+4J 
(a is the lattice spacing) and the operator correspondences ~ & x , ~ ° v ■ We note 
that although continuum limits of H have been studied earlier ||, the identification of the 
universality class of Mi, M 2 has not been made. 

The required correlators can be obtained from a scaling analysis of results of Barouch 
and McCoy |§ on finite T correlators. Express a x in terms of the Jordan- Wigner fermions, 
and evaluate the resulting fermion correlators; this yields an expression in the form of a 
Toeplitz determinant [|7|,|| 



(Tfaf^) = det \Dl- m \ l= i... n>m= i... n (11) 

where D e = ^(d^^e-^D^) with D(<f>) = e^tanh [\E(<f>)\/k B T\ (£(0)/| £(</>) I) and 
E(4>) = 2 J cos + A — i2J / -fsin<f). It now remains to take the large n limit of flTID ; in 
general, this is quite difficult, and leads to a computation of considerable complexity ||. 
However, in a limited portion of the phase diagram (|A/2J| < 1, 7 > which corresponds 
to x > 0, y > in the scaling limit associated with Mi) the computation is simpler because 



it is possible to directly apply Szego's lemma flQ|| . This limited result is all we shall need 



as it is possible to deduce the scaling functions elsewhere by the powerful requirement that 
both F and G are analytic for all real, finite x and y ||11|| . This analyticity is a consequence 
of the absence of thermodynamic singularities in one dimensional quantum systems at any 
finite temperature. The use of analyticity was essential in our being able to express the final 
results in a compact form, and it would have been practically impossible to see the hidden 
structure in the very lengthy results of Ref || otherwise. 
In its region of applicability, Szego's lemma [HJ tells us 



hm (<rf <xf +n ) ~ e™° exp I E W*-*J ( 12 ) 

where ha. D ((f)) = Y^^L-oo oi v e %v ^ . We can read off results for £ and A by comparing (^) with 
(0). We obtained for the scaling function of the correlation length 

t?( \ ^ d< i^ ^ [qV + (x - q 2 ) 2 } 1/2 , al s, ,, af , (y 2 -4x)V 2 - \y\ 

F(x,y)= — lncoth + 9(-y)\y\ + 9(-x) (13) 

Jo n Z Z 



Only the x > 0, y > portion of flT3| ) was obtained from (|T2|); the remainder was deduced 



by the requirement of analyticity. Indeed, even though is appears otherwise, the result ( |13| ) 
is in fact a smooth, differentiable function of x, y for all real x, y including along the lines 
x = and y — 0. 

It is quite interesting to see how the Ising and (Ising) 2 behavior of the critical lines in 
Fig |l] emerges from (|13"D. First, we observe that along the line A 2 = 

f B (x)=F(x,0) (14) 

where fs precisely the correlation length crossover function for the U(l) invariant dilute 
Bose gas in the form presented in |4 , 11] , where it was obtained from the results of Ref [12 



The Ising transition realized by crossing the fi = axis for finite A 2 > is obtained as 
follows 

fi(x) = lim yF(x,y), (15) 

y »oo 

where // is now the correlation length crossover function of the transverse-field Ising model 



obtained in Ref [fLl|| . The prefactor of y on the r.h.s. of fllq) ensures that // is multiplied 



by a power of T appropriate to the z = 1 Ising transition. Crossing fi = axis for A2 < 
requires one to compute (^y^y) to obtain //. Finally, the (Ising) 2 transition realized by 
crossing the A2 axis for fi > is characterized by the limit 

/j(s) = Jim x 1/2 F{x,y = s/x 1 ' 2 ). (16) 

Precisely the same function // emerges in the very distinct limits in QT5D and (0). 

For the amplitude of the correlation function, the computation of G(x > 0,y > 0) from 
(HI) Hi l ea d to a very lengthy and complicated result. However, it was found that the 



analogous result for the Ising model [I1J simplified considerably when expressed in terms 
of derivatives of the correlation length crossover function. We found the same remarkable 
simplification here, and the analytic continuation to all x, y was then straightfoward; we 
obtained 



too , \J dF(x\y) \ 2 1 / dF(x',y) dF(x',y) \ 2 \ f^ 2 dx' 

iog[G{*,v)) = -J x dx^—^-j -^{y^r- + ^-) )-L ^ 

(17) 

Both integrands above are singular at x' = 0, but the singularities cancel in the sum; 
it can be verified that the expression (|T7|) defines G as function which is smooth for all 
real x^y ^ cis required. The limiting behavior of G near the critical lines of Fig [l] is sim- 
ilar to that of F\ we have Qb{x) = G(x,0) which is the analog of ( |14D (with gg now 
the crossover function of the amplitude of the U(l) invariant dilute Bose gas [[TT| ,|r2|j ) , 
log[(7/(x)] = lim^oo \og{(2/yy/ 2 G(x, y)] as the analog of ( |T5| ) for Ising transition line (with 
gi the crossover functions of the amplitude of the Ising model 0]), and log[^/(s)] = 
(1/2) lim^^oo log[2 1/,2 G(x, y = s/x 1 / 2 )} as the analog of ([IB]) for the (Ising) 2 transition. Es- 
tablishing these limits required use of the recently introduced ]TTJ identities obeyed by 
Glaisher's constant. 

From the above finite T results for the amplitude, we can also deduce the value of the 
spontaneous magnetization: limT-+o^4 = (^x) 2 in a regime where liniT^o^ -1 — 0. This 
method gives us a T = universal function $, 
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<**> = e^)6{\ 2 )(2m^l% f^T^) (18) 

with $(s) = (s/2) 1//4 , which describes the crossover of the spontaneous magnetization be- 
tween the Ising and (Ising) 2 lines of Fig |l|. 

Finally we note that exact, finite temperature, crossover functions of spin correlators 
near bulk quantum critical points, below their upper critical dimension, had previously 
been computed only for the d = 1 dilute Bose gas |T^,||,|TTJ and the d = 1 transverse field 
Ising model JTTJ ; our results flT3| ) and QT7| ) for F, G contain these as earlier results as limiting 
cases, and also (universally) interpolate between them. 
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FIG. 1. Phase diagram of Sf (equivalent to Sq + £2 in d = 1) at T = 0. The lines denote 
second-order quantum phase transitions with gapless excitations, and there is gap to all excitations 
elsewhere. The multicritical point M is at fi = 0, A2 = and has z = 2. The action has U(l) 
invariance only for A2 = 0, and the conserved U(l) charge, (^ b ^>b) = (^f^f), is quantized at 
zero for A2 = 0, \i < 0. On the critical line \i > 0, A2 = 0, this charge varies continuously and the 
system is described at sufficiently long scales by the conformally and U{\) invariant (Ising) 2 theory 
with z = 1, central charge c = 1. Similarly, the critical lines \x = 0, A2 > and [i = 0, A2 < 
eventually map on to the conformally invariant {z = 1) Ising field theory with c = 1/2. The point 
M is the only truly scale invariant point, but it is not conformally (or even Lorentz) invariant. 
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FIG. 2. Phase diagram of the lattice model if at T = 0. The expectation values of a x and a y 
vanish unless otherwise noted. The point M\ is at 7 = 0, A/J = —2, while the point M2 is at 
7 = 0, A/J = 2. The vicinities of both Mi and M% map independently onto the continuum model 
of Figure [j]. 
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